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1. Introduction 

The study of QCD at non-zero baryon density by numerical simulations on a space-time lattice 
is plagued by the well-known sign problem: the fermion determinant is complex and the Monte 
Carlo sampling becomes unfeasible. One of the possibilities to circumvent this problem is to 
perform Monte Carlo numerical simulations for imaginary values of the baryon chemical potential, 
where the fermion determinant is real and the sign problem is absent, and to infer the behavior at 
real chemical potential by analytic continuation. The idea of formulating a theory at imaginary /i 
was first suggested in Ref. while the effectiveness of the method of analytic continuation was 
pushed forward in Ref. [^. Since then, the method has been extensively applied to QCD [^- IC] 



and tested in QCD-like theories free of the sign problem [ 11 -[T^ and in spin models [17, 18]. The 
state-of-the-art is the following: 

- the method is well-founded and works fine within the limitations posed by the presence of 
non-analyticities and by the periodicity of the theory with imaginary chemical potential [ [T9| ] ; 

- the analytic continuation of physical observables is improved if ratios of polynomials (or Pade 
approximants [20]) are used as interpolating functions at imaginary chemical potential [|l^, 14]; 



- the analytic continuation of the (pseudo-)critical line on the temperature - chemical poten- 
tial plane is well-justified, but a careful test in two-color QCD [ p^ ] has cast some doubts on its 
reliability. 

In particular, the numerical analysis in two-color QCD of Ref. [|l^] has shown that, while there 
is no doubt that an analytic function exists which interpolates numerical data for the pseudo-critical 
couplings for both imaginary and real /i across /X = 0, determining this function by an interpolation 
of data at imaginary /i could be misleading. Indeed, in the case of polynomial interpolations, there 
is a clear indication in two-color QCD that non-linear terms in jx^ play a relevant role at real jX, but 
are less visible at imaginary /i, thus calling for an accurate knowledge of the critical line there and, 
consequently, for very precise numerical data. The above described scenario could well be peculiar 



to two-color QCD and strongly depend on the choice of parameters of Ref. [14]. Therefore, in 
this work we perform a systematic study of the analytic continuation of the critical line in another 
sign-free theory, SU(3) with a non-zero density of isospin. The dependence on the fermion mass 



in two-color QCD is considered in a separate contribution [piQ. The aim of this study is to single 
out some general features of the analytic continuation of the critical line and to understand if and 
to what extent they can apply also to the physically relevant case of QCD. 

2. Analytic continuation of the critical line in three-color QCD at finite isospin 
chemical potential 



Three-color QCD with a finite density of isospin charge [ [22| is a theory in which the chemical 
potential is pL for half of the fermion species and —pL for the other half. The partition function, 
which is even in pi and depends only on pL^, can be written as follows: 

Z{T,ix) = J me-^°detM[n]detM[-^] , (2.1) 

where the integration is over gauge link variables. So is the pure gauge action and M the fermion 
matrix (we adopt a standard staggered discretization). This leads to a real and positive measure. 
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Figure 1: (Left) Distribution of the real part of the Polyakov loop in SU(3) with finite isospin density on a 
8^ X 4 lattice with am=Q.l at /x^/ (;rr)^=0.16 and for three )3 values around the transition. (Right) Critical 
couplings obtained in SU(3) with finite isospin density on a 8^ x 4 lattice with am=0.\, together with a 
polynomial fit of order ii^ to all data. 



because of the property detM[— /x] = (detM[/i])*, and therefore to a theory free of the sign problem. 
This theory is obviously closer to real QCD than two-color QCD, being yet unphysical, since it 
implies a zero baryon density, while in Nature a non-zero isospin density is always accompanied 
by a non-zero baryon density; moreover the isospin charge is not a conserved number in the real 
world. Nevertheless, for our purposes this theory is very convenient since it provides us with 
another theoretical laboratory for the method of analytic continuation. 

Similarly to SU(2) with finite baryon density, at imaginary values of the chemical potential /i 
the theory exhibits RW-like transition lines, the first RW sector being given by the strip —(0.5)^ < 
pL^ / (tiT)^ < (we refer to Ref. [23] for a detailed discussion of the QCD phase diagram in presence 
of an imaginary isospin chemical potential). 

In our numerical analysis, we consider finite isospin SU(3) with Nf = S degenerate staggered 
fermions of mass am = 0.1 on a 8^ x 4 lattice. The critical line in the temperature - chemical 
potential plane is a line of (strong) first order transitions, over all the investigated range of jj.^ 
values, -0.2304 < lx'^/{nTf < 0.2025. This is one of the reasons for working on a small volume 
(tunneling between the different phases would have been sampled with much more difficulty on 
a larger volume) and clearly emerges from the distribution on the thermal equilibrium ensemble 
of the values of observables like the (real part of) the Polyakov loop, the chiral condensate, the 
plaquette across the transition (see, for example. Fig. [T](left)). Typical statistics have been around 
lOK trajectories of 1 MD unit for each run, growing up to lOOK trajectories for 2-3 j8 values 
around j3c(/i^), for each }JL^, in order to correctly sample the critical behavior at the transition. The 
critical j3 (/x^) is determined as the point where the two peaks have equal height and, in all the cases 
considered, this point turned out to be the same for all the adopted observables. 

2.1 Results for the critical line at finite isospin 

The general strategy is the following: after determining, for a set of }j} values, the critical 
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couplings j3c(/i^), the critical line is guessed by interpolating the values of j8c(/x^) for /x^ < only. 
The validity of the interpolation is evaluated by comparing its analytic continuation to the region 
>0 with the direct determinations of the critical coupling in this region ^ 

We observe from the very beginning that data for j8c(/x^) for both jJ.^ <0 and /x^ > can 
be globally fitted by an analytic function (a polynomial of third order in jj.'^f {hTY nicely works). 
Fig. [I|(right) shows how the fit compares with data. The question is if there are interpolations of the 
critical couplings at /i^ < only, that, when continued to pL^ > 0, agree with the critical couplings 
directly determined in the latter region. 

We have tried several kind of interpolations of the critical couplings at jj.^ < 0. At first, we 
have considered interpolations with polynomials up to order n We can see that data at /x^ < are 
precise enough to be sensitive to terms beyond the order jj.^; indeed, a good ;^^/d.o.f. is not achieved 
before including terms up to the order jj.^, in agreement with the outcome of the global fit discussed 
above. The extrapolation to /i^ > for the polynomial of order /i^ is shown in Fig. §(left); it agrees 
with direct determinations of j8f(/i^), within the 95% CL band. 

Then, we have considered interpolations with ratios of polynomials of order up to n^. In all 
but one cases we got good fits to the data at /x^ < 0, but only two extrapolations to /x^ > compare 
well with numerical data in that region: the ratio of a 4th to 6th order polynomial and the ratio of 
a 6th to 4th order polynomial [23]. It is interesting to observe that the two interpolations which 
"work" have in common the number of parameters. 

Both kinds of fits considered so far have evidenced that the role of terms of order larger than 
jjL'^ cannot be neglected. Since the data more sensitive to these terms are those farther from jj.^ = 0, 
while data closer to /^^ = should "feel" only the /i^ term in a polynomial interpolation, we per- 
formed a fit with a polynomial of the form ao + aipi^ / {^T)^ in a small region jU"^ < and fixed the 
value of the parameter ai (see Ref. p3[] for details). Then, we kept ai fixed and repeated the fit 



'We refer to Ref. [gBJ for all the determinations of the critical couplings in the finite isospin SU(3) theory on a 
8"' X 4 lattice with fermion mass am=0. 1 and for the parameters of all the fits presented below. 
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on all available data at /x^ < with a polynomial of the form + a\}JL^ / (tiT)'^ + a2/i^/ (tiT)^ + 
a-ijX^ / {nT)^ . The resulting interpolation and its extrapolation to the region /x^ > are shown in 
Fig. ^(right). The comparison with critical couplings at /i^ > is good and the 95% CL band is nar- 
rower than in the unconstrained /x^-polynomial fit (see Fig. fright)), meaning that this procedure 
leads to increased predictivity for the method of analytic continuation. 

At last, we have attempted the fit strategy to write the interpolating function in physical units 
and to deduce from it the functional dependence of j3c on /i^, after establishing a suitable corre- 
spondence between physical and lattice units. The natural, dimensionless variables of our theory 
are T /Tc{0), where Tc{0) is the critical temperature at zero chemical potential, and ijl/{hT). The 
question that we want to answer is if fitting directly the dependence of T/Tc{0) on }x/{kT) may 
lead to increased predictivity for analytic continuation. We shall name this kind of fits as "physical" 
fits. While ijl/{hT) is one of the dimensionless variables used in our simulations, T/Tc{0) is not 
and must be deduced from the relation T = 1/ {Nta{l5)), where Nt is the number of lattice sites in 
the temporal direction and a{p) is the lattice spacing at a given j3 ^. Since our determinations for 
range between ~ 4.5585 and ~ 4.842, it can make sense to use for ^(jS) the perturbative 2-loop 
expression with A'^^. = 3 and Nf = S. 

We have tried several different fitting functions and report two cases which work particularly 
well. The first is given by the following 3-parameter function: 

'Tcitl 
Tc{0) 

^Strictly speaking the lattice spacing depends also on thie bare quark mass, whichi in our runs sliglitly ciianges as 
we ciiange jS since we fix am. However in the following evaluation, which is only based on the perturbative 2-loop 
/3 -function, we shall neglect such dependence. 



\+Aii^/{KTAii)Y 



(2.2) 
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(2.3) 



leading to the following implicit relation between jSc and pL^: 

2-loop 2-loop l+-^AtV(?r71'(At))^ 

In Fig. |3|(left) we compare to data the fit with the function (2.3): one can see that the extrapolation 
to the region /x^ > behaves very well. The values of the fit parameters are j8c(0) = 4.6977(13), 
A = -3.25(26) and B = -2.62(12), with ;c2/d.o.f.=1.33. 

As an alternative function for the shape of the critical line, we have tried also the following 

(2.4) 



Tc{0) 




-A)[cos(^)]' 
-A)[cosh(^)]« 



At^ <0 



which explicitly encodes the expected periodicity of the partition function for imaginary /x. The fit 
to data at imaginary /i is very good and its extrapolation to the real chemical potential side compares 
impressively well with data (see Fig.|3|)(right). The resulting fit parameters are j3c(0) = 4.6969(12), 
A = 1.508(15) and B = 0.560(32), with ;t;2/d.o.f.=0.39. This function is a good candidate to pa- 
rameterize the critical line for small values of /i/T. 



In both cases, Eq. (2.3) and Eq.(2.4), the "physical" fit worked very well and with a reduced 
number of parameters with respect to our previous fits, leading to increased predictivity and con- 
sistency with data at real chemical potentials. One can easily check that the adopted functions are 
not appropriate for a continuation of the critical line down to the T = axis, but this is not the aim 
of our study, since such extrapolation would be questionable anyway. 



3. Conclusions 

In this work we have presented results concerning the analytic continuation of the critical line 
in QCD with a finite density of isospin charge. We have detected some features and developed 
some strategies, which could apply and be useful for real QCD at finite baryon density. Let us 
briefly summarize them. 

- Non-linear terms in the dependence of the pseudocritical coupling jS^ on /^^ in general cannot 
be neglected. A polynomial of order jj.^ seems to be sufficient in all explored cases. 

- The coefficients of the linear and non-linear terms in /i^ in a Taylor expansion of j8c(/x^) 
are all negative. That often implies subtle cancellations of non-linear terms at imaginary chemical 
potentials {jx^ < 0) in the region available for analytic continuation (first RW sector). The detection 
of such terms, from simulations at jU^ < only, may be difficult and requires an extremely high 
accuracy. As a matter of fact, the simple use of a sixth order polynomial to fit data at imaginary jx 
leads to poor predictivity, which is slightly improved if ratio of polynomials are used instead. 

- An increased predictivity is achieved if the linear term in /i^ is fixed from data at small values 
of {X^ only. 

- We have proposed a new, alternative ansatz to parameterize the critical line directly in phys- 
ical units in the r,/x plane (instead than in the j8,/i plane) and given two explicit realizations. This 
"physical" ansatz provides a very good description of the critical line, moreover with a reduced 
number of parameters, and leads to an increased predictivity, comparable to that achieved by the 
"constrained" fit. 
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